Evidence for rj c — > 77 and Measurement of J/V 7 — > 37 
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The decay of J/tp to three photons is studied using -i/>(3686) — s- -k + tt~ J/ip in a sample of 106 million 
t/>(3686) events collected with the BESIII detector. First evidence of the decay n c — > 77 is reported, 
and the product branching fraction is determined to be Bj/^,^, 1Vc ^ rll .^ n = (4.5 ± 1.2 ± 0.6) x 10~ 6 , 
where the first error is statistical and the second systematic. The branching ratio for the direct 
decay is = (11.3 ± 1.8 ± 2.0) x 10" 6 . 



PACS numbers: 14.40. Pq, 13.20.Gd, 12.38.Aw 

Decays of positronium to multi-photons are regarded 
as an ideal test-bed for QED Q. The analogous pro- 
cesses in charmonia act as a probe of the strong in- 
teraction 0. For example, the decays of r] c — > 77 
and J/ip — > 37 have relatively simple theoretical de- 
scriptions, and their measurements allow for fundamen- 
tal tests of non-perturbative QCD theory One key 
parameter, which influences the predictive capability of 
QCD models, is the value of the running coupling con- 
stant [H. In principle, the ratios of Bj/^^^/Bj/^^ef, 
and B Vc ->~ ( ~J Bj/^^ee, where Sj/^_j. ee has been precisely 
measured [y] , provide a clean way to extract the running 
coupling constant Here, Bx denotes the branching 
fraction of the decay X. 

With the exception of the CLEO-c measurement of 
Bj/^3~, = (12 ± 3 ± 2) x 10" 6 0], not much is known 
experimentally about the decays of the J/ip to three pho- 
tons. S I?c _ i . 77 is determined mainly from two-photon fu- 
sion 77 M — s- r/ c @, because direct measurements of the 
decay have low statistics. The latest measurement from 
CLEO-c gave Sj/^^ 7 , )c: , )c ^ 77 = (1.2^11 ±0.3) x 10" 6 or 
an upper limit of 6 x 10 -6 at the 90% confidence level 0j- 
Ignoring the higher-order QCD corrections, the branch- 



ing fraction of J/ip — > "fq c ,ri c — > 77 is predicted to be 
(4.4 ± 1.1) x 1CT 6 (|. 

This letter presents the first evidence for rj c — > 77 
in the three photon decay of the J/ip via ^(3686) — > 
tt + ir~ J I 'ip . In addition, the branching fraction of J/ip — > 
37 is measured with improved precision. The analysis is 
based on a sample of (1.06±0.04) x 10 s V(3686) events [| 
collected in the Beijing Spectrometer (BESIII) operating 
in the Beijing Electron-Positron Collider (BEPCII) 
Using ?/»(3686) — > ir + 7r~J/i[> events for this study rather 
than e + e~ — > J jip — > 37 eliminates background from the 
QED process e + e~ — > 37. 

The upgraded BEPCII is designed with a double- 
ring structure. The BESIII detector [l(| is a cylindri- 
cally symmetric detector with five sub-detector compo- 
nents, that are from inside to out: a main drift chamber 
(MDC), a time-of-flight (TOF) system, an electromag- 
netic calorimeter (EMC), a super-conducting solenoid 
magnet, and a muon chamber. The momentum reso- 
lution for charged tracks reconstructed by the MDC is 
0.5% for transverse momenta of lGeV/c. The energy 
resolution for deposited showers in the EMC is 2.5% for 
1 GeV photons. 
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FIG. 1. Distributions of mass recoiling against n + n~ , 
M(7r + 7r~) rcco ji, in data (points), MC simulation of the pro- 
cess J/ip — > 7(77)17 (solid line), and MC 0(3686) inclusive 
decays (dashed line). The arrows indicate the region to select 
J/0 events. This comparison is based on a selected sample of 
J/0 — > 777 — > 37. 

The BESIII detector is modeled with a Monte Carlo 
(MC) simulation based on GEANT4 [Hd- KKMC 
is used to generate MC samples at any specified energy, 
taking into account the initial state radiation and the 
beam energy spread. The -0(3686) decays are generated 
with the generator EVTGEN [l4[ for known decay modes 
according to their branching fractions as listed by the 
Particle Data Group (PDG) 0] , and with Lundcharm [l5[ 
for unknown "0(3686) decay modes. 

For the selection of -0(3686) — > n + n~J/ip, J/ip — ► 3-f 
candidates, events with only two charged tracks and at 
least three photons are required. Minimum distances of 
charged tracks to the interaction point (IP) are required 
to be within 10 cm in the beam direction and within 1 cm 
in the perpendicular plane. The two charged tracks are 
assumed to be n + Tr~ candidates, and the mass recoiling 
against them in the center of mass system must be within 
[3.091,3.103] GeV/c 2 , as shown in Fig.rjJ 

Photon candidates are chosen from isolated clusters 
in the EMC whose energies are larger than 25MeV in 
the barrel region (| cosf < 0.8) and 50MeV in the end- 
cap regions (0.86 < |cos#| < 0.92). Here, 9 is the po- 
lar angle with respect to the beam direction. To reject 
photons from bremsstrahlung and from interactions with 
material, showers within a cone angle of 5° around the 
momenta of charged tracks are rejected. Reconstructed 
showers due to electronic noise or beam backgrounds are 
suppressed by limiting the timing information within [0, 
700] ns after the event start time. 

The 7r + and 7r~ tracks are fit to a common vertex to 
determine the event IP, and a four constraint kinematic 
fit to the initial four-momentum of the 0(3686) is applied 
for each 7r + 7r _ 777 combination. The combination with 
the smallest fit chisquare, xlc ^ s kept, and xlc < 50 is 
required. 

Distributions of M(77)i g versus M(77) sm are shown 
in Fig. [3J where M (77) ig and M(77) sm are the largest 
and smallest two-photon invariant masses among the 
three combinations, respectively. Events from the back- 
ground processes J/ip — > jir / r] / r]' — > 37 are clearly 
seen within the three vertical bands. These backgrounds 
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FIG. 2. Scatter plots of M(77)i g versus M(77) sm for (c) 
data and for MC simulations of the processes (a) J/0 —¥ 
yir°/v/v' -► 37, (b) J/i> -► 7/(1270) -► 7(77)^(77)^, (d) 
J/0 — s> 37, and (f) J/0 — > "/rj c — > 37. Compared with plot 
(c), plot (e) shows the distribution of data after excluding 
backgrounds from J/0 — > jit /rj/rj' , where the empty bands 
are due to the requirements on any two-photon combination. 
The vertical lines indicate the mass windows to reject 7r°, 77 
and 7]' . 

are suppressed by removing all the events where the in- 
variant mass of any two photons, M(77), lies in the 
mass windows [0.10,0.16] GeV/c 2 , [0.50,0.60] GeV/c 2 , or 
[0.90, 1.00] GeV/c 2 . Background events which lie outside 
these mass windows are estimated based on MC simula- 
tion. 

Events from J/ip — > "fe + e~ also contribute to the back- 
ground, when electron-positron tracks fail to be recon- 
structed in the MDC and the associated EMC clusters 
are thus identified as photon candidates. To reject this 
background, the total number of hits in the MDC within 
an opening angle of five EMC crystals around the centers 
of three photon showers is required to be less than 40. 

J/ip — > 77r°7r° events can pass the selection require- 
ments if the two photons from one of the ir° decays are 
nearly collinear or when one of the 7r°s is very soft. Due 
to the large 77r°7r° branching fraction, these events are 
a major source of background. In order to model this 
background, the rich structure of intermediate resonance 
must be taken into account. Hence, a partial wave anal- 
ysis (PWA) [lj] was performed on the BESIII J/ip data, 
following the procedures described in Ref. [l8j]. Figure [3] 
shows excellent agreement between data and MC simu- 
lation based on the PWA results. Decays of J/ip — > "ffj, 
f.J ~ * 77 are negligible because of their extremely small 
branching fractions [||. 

The xlc vame can be used to separate 37 and 77r°7r° 
final states, and the M(77)i g distribution can be used 
to distinguish J/ip — > 7(77)^ from the direct process 
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FIG. 3. Distributions of invariant mass of 7r°7r° (left) and 
polar angle of 7r° (right) in a control sample of 0(3686) — > 
it + tt~ J ftp, J/tp — > 77r°7r°. Points with error bars are data 
and solid lines MC simulation. 



J/tp — > 37. Therefore, a two-dimensional maximum like- 
lihood fit is performed on the M(77)i g and \\c distri- 
butions to estimate the yields of J/tp — > 3j and J/tp — > 
7(77)^. The projections of the two-dimensional fit re- 
sults are shown in Fig. @] and the numerical results listed 
in Table IH The \ 2 per degree of freedom corresponding to 
the fit is 318/349. For the fit, the shapes of both signal 
and background processes are taken from MC simula- 
tion; the normalization of J/tp — > 7(77)^°/?;/^' is fixed 
to the expected density based on MC simulation; and 
the normalization of J/tp — > jn°ir a is allowed to float. 
Backgrounds of non- J/ip decays are estimated using the 
M(7T + 7r-) recoi i sidebands within [2.994, 3.000] GeV/c 2 
and [3.200,3.206] GeV/c 2 . The statistical significance of 
J/tp ->■ 37 {J/tp -> 7(77)^) is 8.3cr (4.1er), as deter- 
mined by the ratio of the maximum likelihood value and 
the likelihood value for a fit under the null hypothesis. 
When the systematic uncertainties are included, the sig- 
nificance becomes 7.3a (3.7a). The branching fraction is 
calculated using B = -rm — ^ 3h3 where 

JV ^(36S6) XO V(3686)^7r+^-J/^ Xe 

n bs is the observed number of events, -/V^( 3686 ) is the 
number of -0(3686) events @, and e the detection effi- 



ciency. B 



V'(3686)- 



r~l~ 71 



is taken from the PDG |f|. 



Many sources of systematic error are considered. 
Simulation of direct J/tp — > 37 assumes the lowest or- 
der matrix element to be similar to the decay of ortho- 
positronium to three photons fl7j . The systematic uncer- 
tainty in the detection efficiency is estimated by weight- 
ing the efficiencies in the Dalitz-like plot of Fig. [U[d) by 
the experimental data. The maximum relative change is 
15%. 

The invariant mass of r\ c in the J/ijj — > jr] c decay is 
assumed to have a relativistic Brcit-Wigner distribution, 
weighted by a factor of E* 3 multiplied by a damping 

factor e -^ 2/8/32 with (3 = (65.0 ± 2.5) MeV [H. Here 
E* is the energy of the radiated photon in the J/tp rest 
frame. An alternative paramctrization of the damping 
factor used by KEDR [20| , as well as variations of the r\ c 
width in the range 22.7-32.7 MeV, are used to determine 
the systematic uncertainties. 

The systematic uncertainty due to possible differences 
between MC and data is evaluated by performing the 
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FIG. 4. (color online) Projection plots of the two-dimensional 
fit to xlc (left) and M(77)i g (right). Here points with error 
bars are data and the solid line is the fit to data. The dotted- 
dash line shows the contributions of J/tp —¥ 37 (dark red), 
the long-dash line J/ip — ► "/rj c — > 37 (red), and the dotted 
line J/tp — > 77r°7T (blue). The shaded histogram represents 
the backgrounds from J/ip — > 771-° /rj/r}' (green) and the solid 
line non- J/ip decays (yellow). 

TABLE I. The detection efficiency e, the signal yields, the 
estimated significance, and the measured branching fractions 
with the errors for the two decay modes. The first errors are 
statistical and the second systematic. Values of the signifi- 
cance outside the parenthesis are statistical only and those 
within the parenthesis include systematic effects. 



modes 


J/tp -> 37 


J/ip -> 777c -> 37 


e(%) 


27.9 ± 0.1 


20.7 ±0.2 


yields 


ff3.4± 18.1 


33.2 ±9.5 


significance 


8.3(7.3)cr 


4.1(3.7)ct 


B(xl0- 6 ) 


11.3 ± 1.8 ± 2.0 


4.5 ±1.2 ±0.6 



two-dimensional fit of the %4C ano - -^(77)ig distribu- 
tions with the MC shapes smeared with an asymmetric 
Gaussian function, the parameters of which are deter- 
mined by comparing a control data sample of J/tp — > jrj, 
■q — ). /y/y to MC simulation. This function serves to adjust 
the detector resolution in the MC simulation to data. 

The efficiency of the mass window requirement around 
the J/tp resonance in the -M(7r + 7r _ ) roco ii spectrum, which 
is sensitive to the tracking efficiency of low momentum 
pions, is studied with a J/tp — > 777, r\ — > 77 control sam- 
ple. The systematic uncertainty in the expected number 
of background events from J/tp —5- jir°(r], rj') was evalu- 
ated by varying their branching fractions by one standard 
deviation [6j. 

It has been verified that the xlc distribution of the 
final states does not depend on the components 
of the intermediate processes involved, i.e., mainly the 
fj states @- However, the M(77)i g distribution does, 
and good understanding of the primary components via 
PWA is needed. Information about the amplitudes in 
J/tp — > 77r°7r° from the previous BESII analysis [l6[ is 
also used in the simulation as an additional check; the 
relative changes in the results are taken as systematic 
uncertainties. 

The photon detection efficiency is studied with dif- 
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FIG. 5. The energy spectrum in the J/ip rest frame of in- 
clusive photons in J/ip — > 37. The points with error bars 
are the partial branching fractions as a function of the ratio 
2E 1 /Mj/^ measured in data. Here E-, is the photon energy 
and Mj/^p is the J/ip mass. The solid line is the theoretical 
calculation according to the ortho-positronium decay formu- 
la H3|. 

TABLE II. Summary of the relative systematic uncertainties. 
A dash (-) means the uncertainty is negligible. 

uncertainties (%) 



sources 


Dj/ip->3i 


^ J j tjj — )■ <y j) c — >■ 3 7 


signal model 


15 




r\c width 




5 


rjc line shape 


1 


1 


resolution 


3 


9 


M (7r + 7r~) rcC oii window 


4 


4 


7T°, rj, rj rejection 


0.5 


5 


PWA model 


2 


2 


photon detection 


3 


3 


tracking 


2 


2 


number of good photons 


0.5 


0.5 


kinematic fit and xlc requirement 


2 


2 


fitting 


5 


5 


number of -0(3686) 


4 


4 


"•0f3686)->7r+ir- J/i/i 


1.2 


1.2 


total 


18 


15 



ferent control samples, such as radiative Bhabha and 
^(3686) TT+TT-J/i;, J/ip -> p°7r° The MDC 

tracking efficiency is studied using selected samples 
of J/ip ->■ pir and 0(3686) -> n + ir-J/i>, J/ip ->■ 
tt + tt~pp [l~8j]- Samples of J/-0 — > 777, ?y — > 77 are selected 
to study the uncertainties arising from requirements on 
the number of photon candidates and the xlc require- 
ment. The uncertainty due to the fitting is evaluated by 
changing the fitting range and the bin width. 



The energy spectrum of inclusive photons in J/ip — > 3j 
provides information on the internal structure of the 
J/ip @ ■ The inclusive photon is defined as any one of the 
three photons in the final state. Partial branching frac- 
tions were measured as a function of inclusive photon en- 
ergy E 1 in the J/ip rest frame. This model-independent 
distribution is illustrated in Fig. [5J where all three pho- 
tons are plotted and the error bars are combinations of 
the statistical and systematic uncertainties. The distri- 
bution agrees well with the theoretical calculation adapt- 
ed from the ortho-positronium decay model. 

In conclusion, the J/ip decays to three photons are 
studied using t/j(3686) — > ir + ir~J/ip decays at BESIII. 
First evidence of the r\ c — > 77 decay is reported, and the 
product branching fraction of J/ip — > 777c and r\ c — > 77 
is determined to be Bj /^^ >lnc ^ c ^ 11 — (4.5 ± 1.2 ± 0.6) x 
1CP 6 . This result is consistent with the theoretical pre- 
diction Q and CLEO-c's result Q- When combined with 
the input of Bj/^^ 7ria = (1.7 ± 0.4) x 10 -2 from the 
PDG 0, we obtain 2?° c ^ 77 = (2.6 ± 0.7 ± 0.7) x 10" 4 , 
which agrees with the result from two- photon fusion . 
The direct decay of J/ip 37 is measured to be 
#7/^37 = (11-3 ± 1.8 ± 2.0) x 10~ 6 , which is con- 
sistent with CLEO-c's result. Considering the ratio of 
Bj/^^zj/Bj/^^ee at the first order of QCD correction, 
the running coupling constant at the J/ip nominal mass 
scale can be extracted. The energy spectrum of inclusive 
photon is also measured. 
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